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IN THE ABSTRACT 

Title: METHOD OF CALIBRATING A SPECTROSCOPIC DEVICE 

The present invention provides a method of calibrating a spectroscopic device for 
providing a non-invasive measurement of an analyte level in a sample. The method comprises 
the steps of: (a) providing a plurality of calibration algorithms; (b) taking a set of non-invasive 
measurements on said sample with said spectroscopic device; (c) calculating a predicted set of 
analyte levels for each of the calibration algorithms in response to the set of non-invasive 
measurements, each of the predicted sets of analyte levels being characterized by a variability 
range, a slope, an ^c^l (a square of the correlation between said set of non-invasive 
measurements and said predicted set of analj^e levels), and a standard error of prediction; and (d) 
selecting an appropriate calibration algorithm by using a suitability score based on the variability 
range, the slope, the R^ and the standard error of prediction for each of the predicted sets of 
analyte levels. A method of generating suitable calibration algorithms in step (a) is also 
provided. 
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METHOD OF CALIBRATING A SPECTROSCOPIC DEVICE 

FIELD OF THE INVENTION 

The present invention relates to the fields of spectroscopy, 
5 spectrophotometry, and chemometrics. In particular, the present 

invention relates to a method of calibrating a spectroscopic device for 
predicting analyte levels in a non-invasive manner. The method of the 
present invention is partictilarly suitable for blood glucose prediction 
based on near-infrared spectrophotometry measurements. 

10 BACKGROUND OF THE INVENTION 

Biotechnological analysis and examination are often based on the 
measurement of various chemical analytes in the composition of a 
biological matrix such as blood, interstitial fluid, or living tissue. Such 
measurements may be used to evaluate a patient's state of health and to 

15 determine what, if any, treatment is necessary. For example, the 

frequent monitoring of blood glucose levels in diabetic persons with 
glucose meters is often necessary to allow such persons to manage the 
diabetes mellitus disease, by taking insulin injections or oral dmgs to 
lower blood glucose when required. Intensive treatment based on 

20 frequent blood glucose measurements can significantly reduce the 

incidence of blindness, kidney loss, and other complications associated 
with diabetes. 

Most home-based glucose measurement systems reqviire the 
patient to invasively collect a blood sample, by pricking his or her finger, 

25 placing the sample on an appropriate test strip, and then testing the 

sample in an optical glucose meter. For millions of diabetics aroiand the 
world, the use of lancets or other sharp instruments to draw blood for 
monitoring tiieir insulin levels is a painful process, and one that often 
builds up calluses on fingers, making tiie collection of blood even more 

30 difficult This invasive procedure may be especially difficult to perform 

on children and therefore partioolarly trying on parents. Furthermore, 
the test strips required for each blood sample are generally not reusable, 
and when multiple measurements are taken each day, amoimt to 
significant costs from the patient's point of view. Thus, despite the fact 

35 that a large nimiber of diabetics should take several measurements 

throughout each day (for some individuals, physicians recommend 
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testing glucose levels from 4 to 7 times daily), due to the pain, cost, and 
inconvenience involved, many diabetics do not monitor their glucose 
levels frequently enough. 

Spectroscopy is based on the analysis of how incident radiation 
5 interacts with the vibrational and rotational states of molecules which 

are of analytical interest. Spectroscopic measurement techniques have 
gained increased popularity becatise of the ability to provide fast zind 
non-invasive measurements of concentrations of different chemicals or 
analytes. For the reasons indicated above, this is particularly desirable 

10 for home based glucose meters. Spectrophotometry is a tjrpe of 

spectroscopy commonly used to quantitatively measxire concentrations 
of analytes based on spectral energy distribution in the absorption 
spectrum of a sample solution or medium. In spectrophotometry, the 
energy distribution is typically analyzed within a range of the visible, 

1 5 ultraviolet, infrared, or near-infrared spectra. For example, near-infrared 

radiation (NIR) is electromagnetic radiation having a wavelengtii of 
between about 0.75 and 2.5 micrometers (Le. from 150 to 400 THz), 
Near-infrared spectrophotometry generally uses instruments with 
quartz prisms in monochromators and with lead sulfide photoconductor 

20 cells as detectors to observe absorption bands, and NIR 

spectrophotometry is increasingly being used to measure in vivo 
analytes sudi as glucose, total hemoglobin, albxmiin, cholesterol, and 
ethanoL 

Non-invasive, spectrophotometric measurement of glucose in 
25 human beings is performed by focusing an incident radiation source (or 

sources) on a specific part of the body and detecting the spectral 
distribution of the radiation transmitted therethrough. The absorbance 
of light from the incident radiation is due to the chemical components 
within that body part including water, fat, protein, hemoglobin, melanin, 
30 and glucose. One difficulty with glucose measurement spectral analysis, 

is the spectral overlap between glucose and other chemicals found in 
blood, often in much greater quantities than glucose. In addition, the 
thickness, color, and structure of ttie skin, bones, and blood through 
which the incident radiation passes wiU affect the transmission (or 
35 reflection/ absorption) thereof. Furthermore, the concentration of 

analytes can vary with changes in activity level, diet, hormone 
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fluctuations, and other factors. Glucose concentration measurements are 
also particularly susceptible to variations in physical and chemical 
conditions including temperature, pressure, humidity, and skin 
hydration. As a result, to perform a reliable non-invasive glucose 
5 prediction, NIR spectral measurements should be performed through a 

vascular equilibrated region of the body, and a NIR glucose 
spectrophotometer must be carefully designed so that tiie quality of raw 
spectral information from an NIR glucose meter is high. See generally 
Waynant and Chenault, "Overview of Non-Invasive Optical Glucose 

10 Morutoring Techiuques", IEEE JJEOS Newsletter, vol. 12, no. 2 (April 

1998); and Bxirmeister and Arnold, "Spectroscopic Considerations for 
Noninvasive Blood Glucose Measurements with Near Infrared 
Spectroscopy", IEEE LEOS Newsletter, vol. 12, no. 2 (April 1998). 

Near-infrared glucose measixrements are generally suitable for 

15 tissue deptiis ranging from about 1 mm to 10 cm, and are often 

performed through a patient's finger tip, although otiier areas of the 
body (for example tiie web between two fingers, an ear lobe, or the 
upper lip) can also be used. The Scimple thickness is an important 
e>^erunental parameter because a greater thickness increases the 

20 amount of absorption and thereby lowers the miiiimum limits for 

detection, whereas because less incident radiation successfully traverses 
through a thicker sample (i.e. without being absorbed) effectively 
increasing the spectral noise: see Burmeister and Arnold, supra. 

In addition to the difficulties of obtcdning acciurate spectral 

25 measurements with NIR spectrophotometry, a more significant 

difficulty associated witti such measurements has been the need to 
calibrate such an instrument so that it may be used by various different 
individuals, whose analyte concentrations and variability, for example, 
may differ substantially. Calibration of spectrophotometers (and of 

30 analytical instruments in general) is necessary to ensure the accuracy of 

measurements performed by such devices. 

Two approaches to calibrate and predict analyte (e.g. glucose) 
levels have been used in the prior art. In the first approach, a thorough 
calibration process is performed separately on each individual with 

35 whom the instrument is to be used. This individual calibration process 

requires taking a ntunber of blood samples over a period of time from 
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the individual, and obtaining reference glucose concentration 
measurements from these samples. A number of corresponding non- 
invasive spectroscopic measurements are taken concurrently, and 
calibration regression analysis is then performed to correlate, e.g. 
5 through linear regression analysis, the non-invasive measurements with 

Ihe reference measurements. In this manner, tiie spectroscopic device or 
spectrophotometer is custom calibrated to the specific individual- A 
major disadvantage of this approach, however, is that the calibration 
model is limited to data from that particular individual which may have 

10 an insufficient amoimt of information regarding variation of spectra 

overlapping (interfering with) those of the analyte of interest. It is 
important to be able to characterize such spectra to enable highly 
accurate analyte predictions. To alleviate this problem, the custom 
calibration process may have to be carried out frequently, even on a 

1 5 daily basis. 

An alternative calibration approach requires that a universal or 
general calibration algorithm be developed. This approach is based on 
the notion that by making a large number of calibration measurements 
(as in the first approach) a single calibration model that sxifficiently 

20 accounts for all individual variability of all overlapping or interfering 

analytes can be calculated. Chemometrics, the application of 
mathematical, statistical and formal logic methods to chemistry, is 
generally used to process and compute the spectral intensity data and to 
produce a calibration model: see Small and Amold, 'T)ata Handling 

25 Issues for Near-Infrared Glucose Measurements", IEEE LEOS Newsletter, 

vol. 12, no. 2 (April 1998); and Shaffer, Small, and Amold, "Genetic 
Algorithm-Based Protocol for Coupling Digital Filtering and Partial 
Least-Squares Regression: Application to the Near-Infrared Anal5rsis of 
Glucose in Biological Matrices," Anal. Chem. 68, 2663-2675 (1996). 

30 However, for measurements of glucose, and otiier ancdytes which vary 

distinctly and differently from one individual sample to another, a single 
universal calibration algorithm is often ineffective and can result in 
significant, and sometimes dangeroxis, erroneous predictions of analyte 
concentration levels. 

35 Thus, while NIR measurements provide a non-invasive, fast, 

painless, and convenient technique to monitor glucose levels, correlation 
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and clinical interpretation of spectral measurements to obtain the true 
glucose levels is crucial for proper therapy and disease management. 
Proper calibration of instruments for different patient populations 
(which will vary in etbnidty, age, weight, and so on) is crucial in 
5 obtaining accurate glucose prediction models. Furthermore, careful 

validation and testing of the non-invasive results and the glucose 
prediction equation is needed to determine if the glucose correlation is 
consistent in all clinically important conditions and for all, or at least 
most, t5rpes of patients. 

10 SUMMARY OF THE INVENTION 

In a first aspect the present invention provides a method of 
calibrating a spectroscopic device for providing a non-invasive 
measiirement of an analyte level in a sample, comprising: 

(a) providing.a pliurality of calibration algorithms; 

1 5 (b) taking a set of non-invasive measurements on said sample 

witti Sedd spectroscopic device; 

(c) calculating a predicted set of analyte levels for each of the 
calibration algorithms in response to the set of non-invasive 
meastirements, each of the predicted sets of analyte levels being 

20 characterized by a variability remge, a slope, an R-^ (a square of the 

correlation between said set of non-invasive measurements and said 
predicted set of analyte levels), and a standard error of prediction; and 

(d) selecting an appropriate calibration algorithm by using a 
suitability score based on the variability range, the slope, the R^ and the 

25 standard error of prediction for each of tiie predicted sets of analyte 

levels. 

Preferably, step (d) comprises a method according to claim 1 
wherein, step (d) comprises 

(i) selecting the predicted sets of analyte levels in which the 
30 standard error of prediction is less than an upper error limit, the 

variability range is greater than a lower range limit, and the slope is 
between a fe&t lower slope limit and an upper slope limit said lower and 
upper slope limits defining an acceptable slope range; 

(ii) for each of the predicted sets selected in step (i), calculating 
35 a suitability score in response to the slope, the R^ and the standard error 

of prediction for that predicted set, and selecting the calibration 
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algorithm corresponding to the predicted set having the optimal 
(highest) suitability score as the appropriate calibration algorithm; 

(iii) if no predicted sets 2ire selected in step (i), selecting the 
predicted sets of analyte levels in which the variability range is lower 

5 than the lower range limit and in which the standard error of prediction 

is less than the upper error limit, 

(iv) from each of tfie predicted sets selected in step (iii), selecting 
the calibration algorithm corresponding to the predicted set having the 
lowest standcird error of prediction as the appropriate calibration 

1 0 algorithm; and 

(v) if no predicted sets are selected in step (i), or step (iii), 
determining that no calibration algorithm is appropriate. 

Also preferably, the acceptable slope range may be subdivided 
into a plurality of ubranges corresponding to a plurality of levels, 
1 5 comprising a first level and one or more subsequent levels, and if no 

predicted sets are selected in the first level, then repeating step (i) at each 
subsequent level imtil a predicted set is selected or there are no more 
subsequent levels. 

Preferably, ttie plurality of calibration algorithms provided in step 
20 (a) are first generated by the steps of: 

(i) compiling non-invasive and corresponding reference data 
sets of measurements of analyte levels for each of a number of samples; 

(ii) rejecting^data sets that are not suitable for calibration; 

(iii) combining data sets that are suitable for calibration into a 
25 plurality of groups depending on whether correlations of the combined 

data sets meet predetermined criteria; and 

(iv) generating a calibration algorithm for each of the groups 
of data sets. 

Steps (a), (c), and (d) of the method of tiie present invention may 
30 be performed on a computer associated with the spectroscopic device. 

Furthermore, ti\e sample may be an individual patient, the spectroscopic 
device may be a near-infreured spectrophotometer, and the analyte may 
be selected from the group consisting of glucose, hemoglobin, albixmin, 
cholesterol, and ethanol. 
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The objects and advantages of the present invention will be better 
understood and more readily apparent with reference to the remainder 
of the description in conjunction with the accompanying drawings. 
BRIEF DESCRIPTION OF THE DRAWINGS 
5 In the drawings which illustrate, by way of example, a preferred 

embodiment of the invention: 

Figure 1 is a correlation scatterplot of glucose predictions for a 
number of different patients based on a universal, single calibration 
model; 

1 0 Figure 2 is a correlation scatterplot of glucose predictions for a 

number of patients based on the multiple algorithm calibration model 
according to the present invention; 

Figure 3 is another correlation scatterplot of glucose predictions 
for a nimiber of patients based on the multiple algorithm calibration 

1 5 model according to the present invention; 

Figure 4 is a correlation scatterplot of glucose predictions for a 
typical patient based on a universal, single calibration model; and 

Figure 5 is a correlation scatterplot of glucose predictions for the 
same typical patient based on the muMple algorithm calibration model 

2 0 according to the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

As mentioned, the present invention relates to a method of 
calibrating a spectroscopic device for predicting analyte levels in a non- 
invasive manner. The present invention can be used with a typical NIR 

25 spectrophotometer system having a light soiu-ce which is projected 

through the item to be examined, a sample interface mechanism, a 
spectrometer to separate the light into its component wavelengths, a 
detector, amplification electronics and a computer. By measuring the loss 
(absorption), between the source and the detector and applying 

30 appropriate chemometric (mathematical) techniques, it is possible to 

non-invasively determine the chemicals being examined since different 
chemicals absorb different amoimts of light at different wavelengths. 
Such a spectrophotometric device and method are described in detail in 
United States Patent No. 5,361,758, the contents of which are 

35 incorporated herein by virtue of this reference. While the present 

description relates primarily to glucose measurement, one of the major 
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fields of application for NIR measurement at present, it will be 
understood that the principles of the present invention equivalently 
apply to otiier analytes non-invasively measured using various 
spectroscopic techniques. 
5 To utilize ttie NIR spectrum for glucose measxirement, it is 

necessary to use a spectrometer which has wide djmainic range, a high 
signal to noise ratio, and exhibits low scattering losses. The output from 
the spectrometer is used to generate spectra with high precision both in 
absorbance and wavelength. Significant glucose absorption bands are 

1 0 centered about wavelengths of 1.67, 2.13, 2J27, and 233 micrometers (as 

discussed in Small and Arnold, "Data handling Issues for Near-Infrared 
Glucose Measiirements", supra). In addition to these wavelengtiis there 
are significant glucose absorption bands at 0.97 and 1.12 micrometers 
which wavelengths allow transmission measurement through a greater 

1 5 distance in the tissue. To be able to use NIR to measxire a particular 

compoimd/analyte, chemometric mathematical analysis is applied to the 
measured spectrum. The mathematical analysis techniques are carried 
out by a computer equipped v^th advanced software capable of 
interpreting the resulting complex spectra. 

20 To xmiversally calibrate an NIR spectrophotometer, the 

absorbance associated with the compound or analyte of interest must 
first be measured on a relatively large number of samples. These NIR 
measurements axe then compared to measurements made in a more 
traditional and more accurate invasive manner. From these comparisons 

25 an cdgorithm is developed that characterizes the analyte to be measured. 

The methods used to generate calibration models for NIR 
spectrophotometers are sophisticated. An important criterion is tiiat 
calibration samples be uncorrelated vri&i other chemical species 
concentrations. Satisfaction of these criteria is burdensome because of 

30 the large number of chemical species typically present in a biological 

matrix or growth medium. 

The present invention provides a multi-algorithm calibration 
method for use with a spectroscopic device capable of predicting analyte 
levels in a biological matrix, such as blood glucose levels, for a variety of 

35 different samples with significantly greater accuracy tiian prior art 

universal single calibration models. The calibration method of the 
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invention has two key parts: the generation of a plurality of different 
algorithms and the selection of the algorithm which is most appropriate 
for generating an accurate prediction of analyte concentratioru 

To generate suitable calibration algorithms, the non-invasive and 
5 reference measurement data sets for a number of individuals are first 

compiled. Once this is achieved, correlation techniques are carried out to 
relate non-invasive spectral measurements to reference glucose levels. 
At this stage the suitability of individual data samples for calibration 
purposes must be considered. If some of the non-invasive spectral 

1 0 measurements in a data set resiilt in a saturation of the measurement 

equipment, e.g. at the output of an analog-to-digital converter, then the 
data sample should be disqualified for use in a calibration algorithm. 
Similarly, if the range of variability of the non-invcisive spectral 
measxirements in a data set is too small, the data sample should also be 

1 5 disqualified from calibration. Furthermore, if the correlation of glucose 

in a particular data set is xmdesirably correlated to another variable, such 
as the concentration of another chemical, the sample should again be 
disqualified. As a result, a large numbers of sample sets will not be 
suitable for calibration algorithm development. 

20 The grouping or combining of calibration-suitable data sets (each 

specific to a particular individual) into different algorithms occurs on a 
trial and error basis, by repeatedly combining these data sets and 
performing a new correlation on the combined sets of data. If the 
correlations of the combined sets meet certain criteria, then the 

25 combination is acceptable, if not, a different combination is attempted. 

Once a combination is deemed acceptable, the addition of a further data 
set is attempted and its acceptability is determined in turn. In this 
manner, a combination of data sets for an algorithm grows from 2 data 
sets, to 3 data sets, to 4 data sets, and so on imtil a suffidenfly large 

30 number of data sets is obtained from which a reliable calibration 

algorithm can be developed. 

In combining data sets for a single algorithm, two criteria are key. 
First, any correlation of glucose in the combined data set group with 
other parameters must be minimized. Such vindesirable correlations may 

35 be exhibited with different chemical species concentrations or other time 

var3ang quantities. A suitable quantitative test for this criteria is to 
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ensure that all undesirable correlations of the combined data sets, as 
indicated by the square of the correlation coefficient R^^ (for the 

undesirable correlations), are less than a certain limit, such as R^^ < q.IO. 
In known manner, the correlation coefficient, R^, is a measure of the 

5 tendency of two variables to vary together. If the addition of a data set 

to another data set or to a suitable sub-combination of data sets does not 
satisfy this criteria, the new combination will not be acceptable. 

In addition, the combination of data sets can be optimized by 
ensuring that the data sets exhibit a high degree of correlation with one 

1 0 another, so a second criteria in combining data sets should require tiiat 

tiie data sets in a particular group be sufficiently correlated with one 
another. This can be achieved, for example, by ensuring that the 
Standard Error of Prediction (SEP) is minimal and the square of ttie 
correlation coefficient, R^^ together with the slope relating reference 

15 glucose to predicted glucose, are greater tiian about 0.8 for an existing 

sub-combination of data sets and another data set being potentially 
combined with that sub-combination. 

In an altemate embodiment of the present invention, instead of 
successively adding on to, or building up, groups of data sets for a 

20 calibration algorithm, a large group of calibration suitable data sets can 

be initially taken up at random. In this embodiment, data sets are 
removed from the group, and the SEP and correlation criteria above are 
reevaluated. If the SEP, slope and correlation criteria improve after tiie 
removal of the data set, the data set is permanently removed from the 

25 group, whereas if the criteria do not improve, the data set is reinstated 

back into the group. Using this top-down approach, the removal of data 
sets can cease when the criteria meet acceptable limits. 

The SEP, slope, and correlation calculations for determination of 
the calibration algorithm groupings (and the associated calibration 

30 algorithms) become increaisingly complex with increasing amounts of 

data and are preferably carried out using partial least-squares (PLS) 
regression analysis. The PLS analysis techniques are well known to those 
skQled in the art and provide good linear approximation by removing 
considerable redtmdant information: see generally Small and Arnold, 

35 "Data handling Issues for Near-Infrared Glucose Measurements", supra. 
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Other types of analysis, such as partial component analysis (PCA) or 
artificial neural networks (ANN), can also be used. 

This calibration algorithm generation process is used to generate 
at least two groups of data sets (corresponding to groups of individuals), 
5 each providing a calibration model or algorithm for predicting the level 

of analyte present in a biological matrix based on a non-invasive 
spectroscopic measurement. Preferably, two initial groups of data set are 
developed, and, from tiiese two groups, additional groups can be 
generated. These additional groups are generally smaller in size than the 
1 0 two initial groups and may meet more stringent correlation criteria. 

These additional groups may be constructed with the top-dovm 
approach described above, wherein data sets are selectively removed 
from an initial group to determine if the correlation criteria improve or 
deteriorate after the removal. 
15 Given the availability of multiple calibration algorithms, it is 

necessary to select the most suitable of these algorithms for predicting 
the levels of analyte in any given patient or sample. Once this is done the 
spectroscopic device is calibrated v^th the selected algorithm. The 
algorittim selection process is generally independent of the process used 
20 to generate the multiple algorithms, although it will clearly depend on 

the number of different algorithms that are available. 

For any given patient or sample on which the non-invasive 
spectroscopic measurement is to be made, it is necessary first to compile 
a set of hig^y accurate reference measurements and corresponding 
25 non-invasive spectral measurements at selected intervals over a period 

of time. In a preferred embodiment, two successive non-invasive 
spectral measurements are successively taken each measiurement 
interval, and their mean is then taken as a single non-invasive 
measurement for that interval. 
30 Before algorithm selection begins, two variables may be set to 

determine the stringency of acceptance requirements for patients (i.e. 
whether at least one of the available calibration algorithms will be found 
suitable for a patient). These are (i) a lower range limit for the range of 
glucose or analyte variability in the data set (if the data does not vary 
35 significantly it may not provide sufficient iivformation to be useful), and 

(ii) an upper error Hmit (SEPj^^x) standard error of prediction 



WO 0m6579 PCT/CAOO/01005 



- 12 - 

(SEP), i.e. the square root of the mean of the squared deviations of the 
reference analyte concentrations from the analyte concentration values 
predicted by a calibration algorithm- 

From the non-invasive spectral measxarements, a predicted set of 
analyte concentration levels are calculated according to eadi calibration 
algorithm. When compared to the reference set of concentration levels, 
each predicted set will be calculated from a linear regression or •'best fit" 
line characterized by; a range of glucose level variability; a slope (which 
ideally equals 1, i.e. the predicted values correspond directly to the 
reference values); a correlation coefficient R, defining the correlation of 
the predicted value set with the reference value set; and a standard error 
of prediction (SEP), as defined above. Note that as the correlation 
coefficient R racreases the SEP will decrease, and vice versa, except in 
cases in which there is an offset or bias in the predicted analyte. 

By way of example, and not to restrict the scope of the present 
invention, the algorithm selection process can include 4 levels or steps. In 
a first level, if the predicted data set for an algorithm meets the following 
criteria 

i) SEP < upper error limit (SEPj^^^) 

ii) Analyte range > lower range limit 
and iii) 0.5 < slope (of regression line) < 1.05, 

a suitability score is calculated according to the follovdng equation: 

score = (slope) (R^) / (SEP) 

If more than one algoritihm meets ttie firet level criteria, the algorittun 
with the highest suitability score is selected and the selection process is 
complete. If no algoritiim qualifies in terms of the first level criteria, the 
selection process advances to a second level. 

The second level is identical to the first level, except that ttie slope 
criteria becomes 



0.33 < slope (of regression line) < 0.5 
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For algorithms diat meet this revised slope criterion (and the SEP and 
range criteria), a sviitability score is calculated as in the first level. Again, 
if more than one algorithm qualifies at ttie second level, the algorithm 
with tixe hi^est score is selected, and the selection process is complete. 

At a third level, if none of the algorithms qualify at either the first 
level or the second level and the range of analjrte measurements is 
greater than the lower range limit, the sample or patient is excluded 
from measurement, since imder these circumstances no calibration 
algorithm can be relied on to perform accurately. 

Lastiy, at a fourtti level, if none of the algorithms qualify at eidier 
the first level or the second level and the range of anal3rte measurements 
is lower iiian the lower range Umit, the algorithm with the lowest SEP is 
selected if that SEP is less than the upper error limit If the SEP is greater 
than the upper error limit, the patient is again excluded from 
measurement. Thus, in the fourth level, the SEP becomes the sole factor 
in deciding whether an algorittun should be accepted. 

Note that the slope criteria limits in the first and second levels 
described above are preferred limits, and these levels can also be varied 
somewhat (similar to the lower range limit and upper error limit) 
depending on the application and acceptance requirements. Also, the 
number of levels or steps used in the algorithm selection process can 
vary, from a single level covering the entire preferred range, to two or 
more levels subdividing the preferred range into two or more range 
limits or range criterion. Regardless of the nxmiber of levels, however, 
the third and fovurth levels as described above will correspond to the last 
two levels, in which the patient is either excluded, or a decision to accept 
an algorithm is made on the basis of SEP alone. It will be readily 
appreciated by those skilled in the art that the number of levels or steps 
chosen is a matter of preference only. 

To illustrate the principals of the present invention, the multiple 
algorithm calibration method was used to calibrate a NIR 
spectrophotometer for measuring glucose levels in a number of 
different patients. All reference data measurements were taken with a 
glucose analyzer such as the industry standard YSF" (Yellow Springs 
Instruments) Glucose Analyzer which provides precise glucose 
measurements on (invasive) blood samples. The inventor found that 
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only approximately 6% of single NIR patient data sets qualified for 
calibration purposes, givm the criteria described above. These 
calibration-qualified sets of data were grouped into two initial calibration 
cdgorithm sets, each initial group was built up to a group with 33 sets of 
5 data within it. Eleven additional groups were spun off from these two 

initial groups and each of these groups had at least 8 sets of data within 
it. Thus, thirteen different calibration algorithms were available for 
predicting glucose levels. 

The algorithm selection was based on a total of 24 reference 

1 0 measurements (i.e. YSP*^ measured) and 24 means of two NIR finger 

measvurements, taken from a patient over a diree day period. The above 
described algorithm selection process was rxm for a nimiber of different 
patients (not involved in the calibration process). 

Figure 1 shows a correlation scatterplot of glucose predictions for 

15 a number of different patients based on a universal calibration model 

(QUAD 37). The scatterplot ia Figure 1 is superimposed on a Clarke 
error grid, as described in Clarke and Cox, "Error Grid Analysis", 
Diabetes Care, 10:622-628 (1987). The Clarke error grid breaks the 
correlation space into five regions (A-E) that assess measurement 

20 accuracy on the basis of validity of the corresponding clinical decision. 

Correlation points falling within the "A" region correspond to the correct 
clinical decision being made based on the similarity between the actual 
and predicted glucose levels, (i-e. the predicted values deviate by no 
more than 20% from the reference values). In the *'B" region, the 

25 predicted values deviate by more than 20% from the reference values, 

but treatment decisions made bcised on the predicted levels of glucose 
would not jeopardize or adversely afiect the patient. In regions "C", "D", 
and "E", the predictions significantiy deviate firom the reference values, 
and treatment decisions based on these predictions may well be harmful 

30 to a patient. The Clarke error grid is often used to evaluate tiie cUnical 

consequences of home-use blood glucose monitor errors, in contrast to 

more conventional analytical methods which may be more suitable for 

laboratory reference devices. 

Referring to Figure 1, it can be seen that the upper error limit 
35 S^max *^ universal calibration algorithm is 3.07 mmol/L and that 

significant amounts of the data are outside acceptable limits or regions. 
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In contrast. Figure 2 shows the correlation scatterplot of glucose 

predictions for 25 patients based on the multiple algorithm calibration 
model according to the present invention, with SEPj^^^^ equal to 2.76 

mmol/L. Similarly, Figure 3 shows a correlation scatterplot of glucose 

predictions for 22 patients based on the multiple algorithm calibration 
model according to the present invention, witfi SEPj,^^ equal to 2.28 

nunol/L. Both the predictions in Figure 2 and particularly Figure 3 
exhibit much better accuracy than those in Figure 1 (calibrated according 
to a universal algorithm). 

Similarly, Figure 4 shows a set of correlation data for a typical 
patient, calibrated according to the same universal algorithm, wheresis 
Figure 5 shows the correlation data for the identical patient calibrated 
according to the midtiple algoritiim method of the present invention. 
Once more, the improved accuracy of the glucose predictions based on 
the calibration method of the present invention is evident- 
It will be clear to those skilled in the art that the principles of the 
present invention are applicable to a broad range of spectroscopic 
applications. For example, NIR radiation can be used to penetrate such 
items as human tissue, vials of blood, or containers of milk - all normally 
considered opaque. Mathematiccil analysis of the resulting absorption 
spectrum determines the composition of the substance penetrated by 
the light Thus, in a blood testing lab, samples can be examined without 
contact by passing the lig^t tiirough a plastic vial containing the liquid. 
Similarly, in a dairy application, the NCR light is passed through a bag of 
milk and such parameters as butterfat, solid particles and lactate are 
measured. The multiple algorithm method of the present invention is 
easily extendible to these types of non-invasive measurements as well. 
Furthermore, the method of the present invention can also be applied to 
the non invzisive measurement of a nimiber of other blood analytes such 
as cholesterol, hemoglobin, HfclAc, Fructosamine, and 1.5 AG. 

While preferred embodiments of the present invention have been 
described, the embodiments disclosed are iQustrative and not restrictive, 
and the invention is intended to be defined by the appended claims. 
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I CLAIM: 

1. A method of calibrating a spectroscopic device for providing a 
non-invasive measurement of an analyte level in a sample, comprising: 

5 (a) providing_a plurality of calibration algorithms; 

(b) taking a set of non-invasive measurements on said sample 
v^th said spectroscopic device; 

(c) calculating a predicted set of analyte levels for each of the 
calibration algorithms in response to the set of non-invasive 

10 measurements, each of the predicted sets of analyte levels being 

characterized by a variability range, a slope, an (a square of the 
correlation between said set of non-invasive measurements and said 
predicted set of analyte levels), and a standard error of prediction; and 

(d) selecting an appropriate calibration algorithm by using a 
1 5 suitability score based on the variability range, the slope, the R-^ and the 

standard error of prediction for each of the predicted sets of analyte 
levels. 

2. A method according to claim 1 wherein, step (d) comprises 

20 (i) selecting the predicted sets of analyte levels in which the 

standard error of prediction is less than an upper error limit, the 
variability range is greater than a lower range limit, and the slope is 
between a first lower slope limit and an upper slope limit said lower and 
upper slope limits defining an acceptable slope range; 

25 (ii) for each of the predicted sets selected in step (i), calculating 

a suitability score in response to the slope, the R^ and the standard error 
of prediction for that predicted set, and selecting tiie calibration 
algorithm corresponding to the predicted set having the optimal 
(highest) suitability score as the appropriate calibration algorithm; 

30 (iii) if no predicted sets are selected in step (i), selecting the 

predicted sets of analyte levels in which the variability range is lower 
tiian the lower range limit and in which the standard error of prediction 
is less than the upper error limit, 

(iv) from each of the predicted sets selected in step (iii), selecting 

35 the calibration algoritiim corresponding to the predicted set having the 
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lowest standard error of prediction as the appropriate calibration 
algorithm; and 

(v) if no predicted sets are selected in step (i), or step (iii), 
determiiung that no calibration algorithm is appropriate. 

5 

3. A method according to claim 2 wherein, said acceptable slope 
rzinge is subdivided into a plurality of subranges corresponding to a 
pluraliiy of levels, comprising a first level and one or more subsequent 
levels^ and if no predicted sets are selected in the first level, then 

1 0 repeating step (i) at each subsequent level imtil a predicted set is selected 

or there are no more subsequent levels. 

4. A method according to claim 2 wherein, said lower slope limit is 
less than one and said upper slope limit is greater than one- 

15 

5. A method according to claim 4 wherein, said lower slope limit is 
about 0.3 and said upper slope limit is about 1*05. 



6. A method according to claim 3, wherein, said acceptable slope 
20 range is defined by a lower slope limit of about 0.3 an upper slope limit 

of about 1.05. 

7. A method according to claim 1 wherein, the plvtrality of 
calibration algorithms provided in step (a) are first generated by the 

25 steps of: 

(i) compiling non-invasive and corresponding reference data 
sets of measurements of analyte levels for each of a number of samples; 

(ii) rejecting data sets that are not suitable for calibration; 

(iii) combining data sets that are suitable for calibration into a 
30 plurality of groups depending on whether correlations of the combined 

data sets meet predetermined criteria; and 

(iv) generating a calibration algorithm for each of the groups 
of data sets. 



35 



8. A method according to claiim 7 wherein the predetermined 
criteria in step (iii) include rnirdrnizing correlations of the combined data 
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sets in a particular group with parameters other than said analyte and 
maximizing the correlation between data sets in a particular group. 

9. A method according to claim 7 wherein steps (iii) and (iv) are 
5 performed using partial least-squares regression anal3^is- 

10. A method according to claim 7 wherein steps (a)^ (c), and (d) are 
performed on a computer associated with the spectroscopic device, 

10 11. A method according to any one of claims 1-10 wherein the sample 

is an individual patient, the spectroscopic device is a near-infrared 
spectrophotometer, and tiie analyte is selected from the group 
consisting of glucose, hemoglobin,, albumin^ cholesterol, and ethanoL 
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(57) Abstract: The present invention 
provides a method of calibrating a 
spectroscopic device for providing a 
non-invasive measurement of an analyte 
level in a sample. The mediod comprises 
the steps of: (a) providing a pliuality of 
calibration algorithms; (b) taking a set 
of non-invasive measurements on said 
sample with said spectroscopic device; 

(c) calculating a predicted set of analyte 
levels for each of the calibration algorithms 
in response to the set of non-invasive 
measurements, each of the predicted sets 
of analyte levels being characterized by 
a variability range, a slope, an Ri2i (a 
square of the correlation between said 
set of non-invasive measurements and 
said predicted set of analyte levels), 
and a standard error of prediction; and 

(d) selecting an appropriate calibration 
algorithm by using a suitability score based 
on the variability range, the slope, the 
and the standard error of prediction for 
each of the predicted sets of analyte levels. 
A method of generating suitable calibration 
algorithms in step (a) is also provided. 
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3. was filed on as application serial no, and was amended on (if applicable) (in the case of a PCT-filcd application) 
lelOTbed and claimed in international no. PCT/CAOO/01005 filed August 31, 2000 and as amended on (jUTany), which I have reviewed 
md for which I solicit a United States patent. 

[ hereby state that I have reviewed and understand the contents of the above*identified specification, including the claims, as amended by 
any amendment referred to above, 

( hereby claim foreign priority benefits under Title 35, United States Code, § 1 19/365 of any foreign applicatjon(s) for patent or inventor's 
certificate listed below and have also identified below any foreign application for patent or inventor's certificate having a filing date before 
that of the application on the basis of which priority is claimed: 



a. Q no such appUcations have been tiled. 

b. 1571 such applications have been filed as follows: 



FOREIGN APPHCAT«ON(S), IV ANY, CLAIMING PRIORITY UNDER 35 USC § 


119 


COUNTRY 


APPLICATION NUMBER 


DATE OF FILING 
(day, month, yenr) 


DATE OF KSSUE 
(day, month, year) 










ALL FORE 


IGN APPLICATION(5), IF ANY, FILED BEFORE THE PRIORITY AFPLIC 


:ationcs) 


COUNTRY 


APPLICATION NUMBER 


DATE OF FILING 
(flay, month, year) 


DATE OF ISSUE 
(day, month, year) 











I hereby claim the benefit imder Title 35, United States Code, § 120/365 of any United States and PCT international application(s) listed 
below and, insofar as the subject matter of each of the claims of this application is not disclosed in the prior United States application in the 
manner provided by the first paragraph of Title 35, United States Code, § 1 12, 1 acknowledge the duty to disclose material information as 
defined in Title 37, Code of Federal Regulations, § l-56(a) which occurred between the filing date of the prior application and the national 
or per international filing date of this application. 



V,S. APPLICATION NUMBER 


DATE OF FILING (day, month, year) 


STATUS (patented, pending, abandoned) 









I hereby claim the benefit under Title 35, United States Code § 1 19(e) of any United States provisional application(s) listed below: 



U.S. PROVISIONAL APPLICATION NUMBER 


DATE OF PILING (Day, Month. Year) 


' : 60/151,536 


31 August 1999 
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ackiiowiedge &c Suty to disclose information that i'a material to the patentability of this application in accordance with Title 37, Code of 
ederal JLegulations, § 1-56 (reprinted below): 

i 1.56 Duty to disclose information material to patentability. 

(a) A patent by its very nature is affected with a public interest. The public itxterest is best served, and the most effective 
jatcnt examination occurs when, at the time an application is being examined, the Office is aware of and evaluates the teachings of all 
nfornwtion material to patentability. Each individual associated with the filing and prosecution of a patent application has a duty of candor 
ixxd good faith in dealing with the Office, which includes a duly to disclose to the Office all information known to that individual to be 
nateriai to patentability as defined in this section. The duty to disclose information exists with respect to each pending claim until the 
;laim is canceled or withdrawn from consideration, or the appUcation becomes abandoned. Information material to the patentability of a 
:laim that is canceled or withdrax^^n j&orn consideration need not be submitted if the information is not material to the patentability of any 
qjaimr e mnjjginj; ^^r^^^ r rnnciHpra ri nn in T h ff applTratinn. There is no duty to submit infonnation which is not material to the patentability of 
my existing claim. The duty to disclose all information known to be material to patentability is deemed to be satisfied if all mfomiation 
oiown to be material to patentability of any claim issned in a patent was cited by the Office or submitted to the Office in the matjner 
prescribed by §§ 1.97(b)-(d) and 1.98. However, no patent will be granted on an application in connection with which fraud on the Of&ce 
ivsis practiced or attempted or the duty of disclosure was violated through bad faith or intentional niisconduct. The Office encourages 
ipplicants to carefully examine: 

« 

( 1 ) prior art cited in search reports of a foreign patent office in a counterpart application, and 

j (2) the closest information over which individuals associated with the filing or prosecution of a patent application 

believe any pending claim patentably defmes, to make sure that any material information contained therein is disclosed to the Office. 

(b) Under this section, information is material to patentability when it is not cumulative to information already of record or 
being made of record in the application, and 

(1) It establishes, by itself or in combination with other information, a priwja facie case of uijpateniability of a claim; 

or ! , 

(2) It refutes, or is inconsistent with, a position the applicant takes in: 

(i) Opposing an argument of unpatentability relied on by the Office, or 

(ii) Asserting an argument of patentability. 

A pritna facie case of ujo^jatentability is established when the information compels a conclusion that a claim is unpatentable under the 
preponderance of evidence, burden-of-proof standard, giving each term in the claim its broadest reasonable constniction consistent with the 
specification, and before any consideration is given to evidence which may be submitted in an attentat to establish a contrary conclusion of 
patentability. 

(c) Individuals associated with the filing or prosecution of a patent application within the meaning of this section are: 

(1) Each inventor named jn the application: 

(2) Each attorney or agent who prepares or prosecutes the application; and 

(3) Every other person who is substantively involved in the preparation or prosecution of the application and who is 
associated with, the inventor, with the assignee or with anyone to whom there is an obligation to assign the apphcation. 

(d) Individuals oiher than the attorney, agent or inventor may comply wifli this section by disclosing information to the 
attorney^ agent, or inventor. 

(e) In any continuation-in-part application, the duty under this secnon includes the duty to disclose to the Office all 
information known to the person to be material to patentability, as defined in paragraph (b) of this section, which became available between 
Che filing date of die prior application and the national or PCT international filing date of the continuation-in-patt application. 
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herefc)r-appoinf the following aitomey(s) and/or patent agent(s) to prosecute this application and to trajosact all business in the Patent and 
'rademark Office connected herewith: 



[}, M. Jeffer 
Utcra, Allan G, 
Vnderson, Oiegg I. 
latzh, Bfian H- 
Jeard, John L- 
iems, John M 
Jranch, John W. 
Jtown, Jef&ey C. 
Jnicss, Steven C. 
jujmham, Jonelle 



Reg. No > 4f 359 ^ 
Reg. No . 4Q,274 
Reg. No . 28,828 
Reg. No- 32^960 
Reg. No;2y,612 
Reg. No:3^S£ 
Reg. No! 41,633 
Reg, No jH,6^ 
Reg. No. -34,13j) 
Rffg N n 4] MO 



Jyme, Linda M. 
^tnpbell, Keith 
!^arlson, Alan G. 
^aispers, Philip P. 
Hhfford, John A, 
^ook, Jeffrey 
Disiignauh, Roaald A. 
Daley, Dennis R. 
)aulTon, Julie R. 
DcVries Smith, Katherine M. 
DiPietro, Mark J. 
Doscotch, Matthew A. 
edell, Robert T- 
Epp Ryan, Sandra 
Fitzsimmons, ICarcn A. 
Goff, Jared S. 
Goggin, Matthew J, 
Gblla, Charles E. 
Gorman, Alan G. 
Gould, John D. 
Gregson, Richard 
Gresens^ John J. 
Hamer, Samuel A. 
Hamre, Curtis B, 
Harrison, Kevin C. 
Hemmxgs, Mark 
Hcrtzberg, Bren A. 
Billson, Randall A. 
Holzer, Jr., Richard J. 
Hope, Leonard J. 
Jardine, John S. 
Johns, Nicholas P. 
Johnston, Scott W. 
Kadievitch, Natalie D. 
Kalinsky, Robert A. 
Kettelbergcr, Denise 
Keys, Jeramie J. 
Knearl, Homer L. 
Kowalchyk, Alan W. 
Kowalchyk, Katherine M. 
Lamberty, Michael 
Lacy, Paul E, 
Larson, James A. 



Reg. No..aZ^4 
Reg, No. 46,597 
Reg, No. 25,9^9 
Reg- NQr33!2?7 
Reg. No. 30.247 
Reg. No. 48,649 
Reg. No'15;96S 
Reg. No. 
Reg. No.j;52l4 
Reg. No. 42.1^ 7 
Reg. N o. 28,70 7 
Reg. No. 48>957 
Reg. No-„2IU67 
Reg. No. 32^7 
Reg. No.JEsSa^TO 
Reg. No, Mh^6 
Reg. No..-44tU5 
Reg.No.-S6^ 
Reg. No. 38,472 
Reg, No. JCllX 
Reg. No- 41.804, 
Reg. No.JiSJJU 
Reg. No. 46,754 
Reg. No r29a6? 
Reg. No. 46J59 
Reg- No-,4si^ 
Reg. ^o.JZMQ 
Reg. No. 31.838 
Reg. No . 42.6£^ 
Reg, No, 44,774 _ 
Reg. No.,4S4^ 
Reg. No.li!99S 
Reg. No- 39.721 
Reg, No-JiJil4 

Rcg.No.P;:52d?l 
Reg. No. 33,924 
Reg. No-'^a^ 
Reg. No. 21,197 
Reg. No- JT2» 
Reg. No. li^P^ 
Reg, No, P3l7g0 
Reg. No. 



Leonard, Christopher J. 
Liepa, Mara E. 
Lindquist, Timothy A. 
Lown, Jean A. 
McDonald, Daniel W. 
Mclntyrc, Jr., William F. 
Mueller, Douglas P. 
Nelson, Anna M. 
parsons, Nancy J. 
Pauly, Daniel M. 



Philhps, John B, 
Pino, Mark J. 
Prendergast, Paul 
Pytel, Melissa J, 
Qualey, Terry 
Randall, Joshua N. 
Reich, John C. 
Reiland, Earl D. 
Schmaltz, David O. 
Schuman, Mark D. 
Schumann, Michael D. 
Scull, Timothy B. 
Sebald^ Gregory A. 
Skoog, Mark T. 
Spellman, Steven J. 
Stewart, Alan R. 
Stoll-DeBell, Kirstin L. 
Sullivan, Timothy 
Sumner, John P. 
Swenson, Erik G. 
Tcllekson, David K. 
Trembath, Jon R. 
Tunheim, Marcia A. 
Underhill, Albert L. 
Vandenburgh, J. Derek 
Wahl, John R. 
Weavei', Paul L. 
Welter, Paul A- 
Whipps, Brian 
Whitaker, John E. 
Wier, David D. 
Williams, Douglas J. 
Withers, James D. 
Wong, Bxyan 
Wong, Thomas S. 
Young, ThoTxias 
ZeuJi, Anthony R. 



Reg, No 
Reg. No 
Reg. No 
Reg. No, 
Reg. No 
Reg, No 
Reg. No 
Reg. No 
Reg. No 
R eg. Nq 



41,94J} 

48.428, 
44,9.? 1, 

40364 



Reg. No. 37,206 
Reg. Nn 4^ f^ ^R 
Reg. No. 4£^ 
Reg. No.dl^ 
Reg. No. 7-^,14g_ 
Reg. No. P^^SO^ 
Reg. N o. 3 7^ 70 3 
Reg. No. -2^67 
Reg. No. 39,828_ 
Reg- No. 31,197 
Reg. No . 30,42 2 
Reg. Nor?37r57 
Reg. No. 325^ 
Reg. No. 40,178 
Reg. No. 45^4 
Reg. No. ^7^074. 
Reg. No.^43J^4 
Reg. No. 42^1 
Reg. No. 29 J 14 
Reg. No. 45.147 
Reg. No. 12,3uU 
Reg, No. aM44 
Reg. No. 42.L89 
Reg. No. 27,40^^ 
Reg. No. ^2T7J 
Reg. No. 33,044 
Reg. No. 4^5^ 
Reg. No. 20.890 
Reg. No- 43,2§1 
Reg. No. 42,222 
Reg. No. 487275 
Reg- No- Tr^'re^ 
Reg. No. 

Reg. No. P-50.836 
Reg- No 48, ^TI^ 
Reg. No.r.?JvZa6 
Reg. No,^45.^i_^ 



Reg. No. 40.44T 

I hereby authorize them to act and rely on instructions from and communicate directly with the person/assignce/attomey/firtn/ organization 
who/which first sends/sent this case to them and by whom/which I hereby declare that I have consented after full disclosure 10 be 
represented unless/until I instruct Merchant & Gould P,C. to the contrary. 

I understand that the execution of this document, and the grant of a power of attorney, does not in itself establish an attorney-client 
relationship between the undersigned and the law firm Merchant & Gould P.C, or any of its attorneys. 
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iease direct all Qprrespoadeucc in this case to Merchant & Gould P.C. at the address indicated below: 

Merchant_& jGiluld 
Minrteapolis, MN 55402-0903 



hereby declare that all statements made herein of my own knowledge are true and that all statements made on inforroation and belief are 
jsUeved to be true; and further that these statements were made with the knowledge that willful false statements and the like so made are 
>uni5hable by fine or imprisonment, or both, under Section 1001 of Title 18 of the United States Code and that such wiUfUl false statements 
aay jeopardize the validity of the application or an5' patent issued thereon. 





Full (Mame 
Of Invtoior 


Family Nftinc 
Cad«ll 


First Given Name 

Theodore ^ 


Second Given Name | 


0 


Residence 
& Citizens lUp 


City 


State or Foreign Country 

Ontario 


Country of OtizenAhSp 

Canada 


1 


Msilirtg ^ 
Addre$9 


Address 

13 Elgin Street East 


City 

Conestogo 


State & Zip Code/Couocry 
Ontario NOB INO / Canada 


signature of Investor ZOl: ^ 






23552 

f AT^hrrTRAPEMARK OFFICE 
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